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ABSTRACT 


The goal of this work is to develop an equivalent circuit mode! for a typical 
discontinuity in shielded suspended substrate line, in particular, the step-change in 
width. This circuit model is needed for use in CAD analysis of microwave systems 
using the suspended stripiine transmission mode. 
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I. INTRODUCTION 


The shielded suspended substrate line (SSSL) is a transmission medium useful 
for radar and circuits at microwave and millimeter-wave frequencies. In order to 
use this transmission medium in the construction of microwave circuits and filters, 
it is necessary to have valid circuit models for typical discontinuities such as the 
series gap in line, open-ended stub, and a discontinuous change in width. 

There is a definite need for an accurate full-wave analysis of strip 
transmission line structures. Even though many studies have been performed for 
the discontinuities of conventional microstrip line, almost no discontinuity 
calculations are available for SSSL. By full-wave analysis is implied the process 
of rigorously solving the frequency-dependent electromagnetic (EM) boundaiy' 
value problem with retention of all the field components. 

In this thesis, the discontinuity structure of interest incorporated in a SSSL 
resonator has been calculated by a fomiulation employing full-wave analysis. The 
solution of the problem has been derived using an efficient method. Specifically, 
the derivation of the characteristic equation for resonant frequencies of the 
resonator model is carried out using Galerkin’s technique applied in the spectral or 
Fourier transform domain instead of the space domain. The resonant frequencies 
of the structures of interest are obtained by numerically solving the characteristic 
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equation. The details of the analysis method will appear in Section II of this 
work. 

In Section III, the step discontinuity will be analyzed by placing the 
discontinuity in an open-ended strip line resonator. In the continuous strip case, the 
strip line has no discontinuity step in width. Following this, the equivalent circuit 
of the discontinuity will be deduced from the numerical resonance data for the 
embedding SSSL resonator when the discontinuity is inserted. 
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II. THEORETICAL ANALYSIS 


The microstrip resonator to be analyzed is shown in Fig.l. A strip 
transmission line which has a step in width is placed on the suspended substrate. 
The shielded suspended substrate line is constructed by placing the substrate within 
a channel milled in the sidewalls. Selection of both depth and width of the shield 
enclosure is important. Dimensions must be sufficiently small to avoid propagation 
of waveguide modes within the shield: these modes would actually produce an 
unwanted coupling with the resonance one is trying to detect. On the other hand, 
if the enclosure dimensions are too small, the propagation characteristics of the 
transmission line itself would be significantly modified. 



Figure 1 End and top view of suspended substrate microstrip line 


3 









It is assumed that the thickness of the conducting strip is negligible and that 
all the media and conductors are lossless. For simplicity, the strip is assumed to 
be located symmetrically in the longitudinal direction within the ends of the 
shielding enclosure. 
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A. FOURIER-TRANSFORMED GALERKIN’S METHOD FOR THE 


SHIELDED STRIP RESONATOR 

The analysis of wave propagation on shielded strip transmission systems of 
the type considered here has been carried out by Itoh and Uwano fRef.2]. These 
authors treated the problem using Galerkin’s method of moments in the Fourier- 
Transform domain. Their work is summarized in terms of a set of Green’s 
function equations: 


In Eqs. (la) and (lb), and are the z and x components of the currents on 
the strip conductor, and and E^ are the components of electric field tangent to 
the substrate surface, and the are the Green’s functions. 

The tildes over the factors in Eqs.(la) and (Ib) imply Fourier transformation 
of the respective quantities. In this work, the Fourier transformation is carried out 
in the bounded region interior to the shielding enclosure placed around the line 
segment representing the resonator. Therefore we have the finite Fourier 
transformation: 
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(2a) 


f (ic^, P^) =£ \f{x,y, z) e^'^'^e^^-^^dxdz 

'T 'T 


where and p„ are the discrete transform variables defined by. 


K = 

(n-l/2)7r/a 

for Ej even, odd(in x) modes 

(2b) 

K = 

rm/a 

for Ej odd, -H, even(in x) modes 

(2c) 

Pn, = 

(m-l/2)7r/RL for E^ even, odd(in z) modes 

(2d) 

P„. = 

mre/RL 

for E, odd, even(in z) modes 

(2e) 


where a is width of waveguide, and RL is length of waveguide enclosure. 

The analytical task then consists of assuming suitable coordinate forms of the 
cunent densities and and Fourier-transforming these. Inner products aie then 
formed in Eqs (la) and (lb), in accordance with the standard procedure of the 
method of moments. This procedure gives rise to a set of homogeneous equations 
in the unknown coefficients assumed in and The solution condition for the 
simultaneous equations leads to a determination of the resonance fr equency of the 
strip line resonator of the problem. The unknowns and E,^ can be eliminated 
by applying Galerkin’s method in the spectral domain. The first step is to expand 
the unknowns and in terms of assumed basis functions and with 
unknown coefficients c„, and d„, . 
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The basis functions and must be chosen to be the Fourier transforms 
of space-domain functions Jzn,(x,z) and 3^n,{XyZ) which are physically realistic, and 
which are zero except for the region kl < w, and bl < f. Substituting (3) into 
(1) and taking inner products of the resulting equations with the basis function and 
for different values of i yields the matrix equations, 



(4b) 

__ w _ ^ _ 

^ XZ^2 ^rr'^ xnT ® i-1,2 , . . , I'l 

m-l n;-' 1 


The right-hand sides of (4) are zero by virtue of Parseval’s theorem, because 
the currents Jxi(x) and the field components Ej(x,h 2 ), E,j(x,h^) vanish in 
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complementary regions of x. For example, if the inner product of E, on the right- 
hand side of Eq.(l) with J^jCkn) is taken, one obtains: 






(x) (-x) dx=0 


In the above, J^,(x) is zero outside the strip, and E^(x) is zero on the strip. 
In this way, the final boundary' condition is satisfied. 

Equations (4) will be expressed in matrix form as follows: 


N M 


E 


1-1,2, . 

.,N 

(5a) 






K 

M 




L 


i = l,2, . 


(5b) 

/ 7 -: 






where from the definition of the inner products associated with the Fourier 
transform defined by (2), the matrix elements are:. 


K, 


, 1 ) 


H'l m-l 


zz'^ ziv P rj 


(6a) 
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(6b) 


K 


( 1 , 2 ) 


n~l m-l 


,IZzxJ 


z-(z.l) 

^2in 


<“o>=EE jz 

n-l /n-1 






(6c) 


(<^c^=EL ^~xip.) 

ti-l nwl 

A homogeneous system of equations is thus obtaine(J in terms of the unknown 
coefficients c„, and dn,. In order that c^, and dn, have nontrivial solutions, the 
determinant of the matrix must be zero, and hence the resonant frequency is 
determined for the resonator and discontinuity structures assumed. 

Equations (5) are now solved for the resonant frequency co^, by seeking the 
root of the resulting characteristic equation. The resonance frequency of the 
suspended stripline resonator is derived from the obtained value of (o„[Eq.(8h)]. 

The accuracy of the solution can be systematically improved by increasing the 
number of basis functions (M,N) and by solving larger size matrix equations. 
However, if the first few basis functions are chosen so as to approximate the actual 
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unknown current distribution reasonably well, the necessary size of the matrix can 
be held small for a given accuracy of the solution, resulting in numerical efficiency. 
Hence the choice of basis functions is important from the numerical point of view. 
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B. GREEN S FUNCTION FOR THE MOMENT METHOD 

Equations (la) and (lb) may be expressed in the matrix form:[Ref.2] 




FJ 12.. 


where the transformed Green’s Functions have the values: 




xs ~ ^ zx~ 


ic^p; 


■2 .-Z^) 


.. ■ "Tl fT^ t-^P 


z. = 


_ y y^^c.^y _ 

’V c. r. -1^+a c. 

r, .. o, -, Y,, Yj,- 
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Y zz^c.*y zi^t- 


y ziy z^C,C,^^y ,_,Y z3C,C,^^y ^,y ,,C,C,^^y‘ 


(6e) 


where the symbols C,.have the values specified below: 
case > 0, 


C,^=coth YjCf, C,^=cothY 2 t, C^^cothy^h (8f) 


case Y,‘ < 0. 


c. 


= -jcoty^d, cot y^t, C^^=-jcoty,h (8g) 


where 




(8h) 


yj 


jwe_j 


( 81 ) 


Y zi 


y X 


(8j) 


where subscripts i refer to the i"* medium. 

The quantities Z„, and Z„ are actually the Fourier tninsfomis of 

dyadic Green’s function components. 
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C. CHOICE OF RESONATOR FREQUENCY BELOW THE CUTOFF OF 
THE SHIELD ENCLOSURE 

The cutoff frequency of the dielectric-loaded waveguide formed by the 
shielded SSSL shown in Fig.l is found to be 5.73GH2, as calculated by the 
equation below [Ref.5]. 


f - -£L t (gj-l) (9a) 

^ ■ 2a^j b e, 

where c is the velocity of light and the value of e, is 2.35. The dimensions a,b 
and t are defined in Fig. 1, and numerical values are shown in Table I. 

TABLE I CONFIGURATION DATA 


Dimension 

Data 

Dimension 

Data 

a 

22.86 mm 

t 

1.59 mm 

b 

10.16 mm 

h 

4.285mm 


5.00 mm 

Er 

2.35 

d 

1 

1 

4.285mm 

RL 

100.00 mm 


Dimensions in Table I correspond to Fig.l as follows: 
a = Width of shield 
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b = Height of shield 

d = Height of lower air layer 

t - Thickr.ess of dielectric substrate layer 

h = Height of upper air layer 

w, = Microstrip line width in continuity case 

e, = Dielectric constant of substrate layer 

RL = Length of shield 

The effective dielectric constant and the frequency determine wavelength 
in the SSSL: 


X 


s 



(9b) 


where c is the velocity of light, f is the frequency, and e^ff is the effective dielectric 
constant. Thus, X^2 determines the fundamental resonance of the open ended 
length of strip, C. 

For simplicity, the strip is assumed to be symmetrically located, although the 
present method of analysis can be easily extended to more general cases. Also, 
the operating frequency is taken to be below cutoff of the shielding waveguide 
partially filled with substrate material, in order to avoid interaction of the strip 
transmission-line resonance with standing-wave resonances of the waveguide modes 
within the shielding enclosure. 
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D. EQUIVALENT NETWORK FOR A DISCONTINUITY IN A 


TRANSMISSION LINE SYSTEM: THE RESONATOR METHOD 

In the resonator method for evaluating an equivalent circuit representation for 
a discontinuity in a transmission-line system, the discontinuity is placed at the 
approximate center of an open-ended strip of the transmission line. 

Ignoring the end effect, the fundamental resonance of the unperturbed 
resonator strip occurs at a total resonator length of XJ2 of the guided wave length 
of the system. For this mode there is an current maximum at the center of the 
resonator (Fig. 2a]. The next higher resonance of the open-ended resonator occurs 
at a resonator length of of the system wavelength [Fig. 2b]. 

In any event, if Che line resonator is housed within a metallic shield, for 
example of waveguide cross-section, it is essential to choose a resonator length 
such that its resonance frequency is well below cutoff of the propagating modes of 
the waveguide system involved. This is a dielectric-loaded waveguide because of 
the presence of the dielectric substrate. 

In the resonator method for determining the equivalent ciicuit representation 
of a discontinuity, the discontinuity is assumed to be lossless, and the equivalent 
circuit is regarded as a two-port circuit element acting as an impedance transformer 
between the two line-segments which make up the resonator[Fig. 3]. 
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Figure 2 Resonance modes 


«l ^ «2 *1 
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Figure 3 Impedances for line segments 


Considering the resonator transmission line segment at the left of the two- 
port, the occurrence of resonance is specified by the condition: 

h = Zinl (10) 

where the star means complex conjugate, z, is the impedance at the input of the 
open stub on the left, and Zj^, is the impedance seen looking into the uiput at port 
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1. The two-port is terminated in the second open-ended line segment which 
makes up the remainder of the resonator[Fig. 4]. 


,_ ^inl 

L... -. 


_1 

1 H 


—-c 


Figure 4 Equivalent-circuit expression for discontinuity 


If the discontinuity under test is physically symmetrical about its center point, 
then it may be represented by a symmetrical T- or PI- circuit[Fig. 5). However, 
for an asymmetrical discontinuity, such as for example the step ci.ange in width, 
the discontinuity may be represented by a general unsymmetrical T or PI, or an L- 
section if appropriate. 

In the evaluation of the discontinuity by the resonance method, a pair of stub 
lengths {| and Cj is chosen, of lengths appropriate to the test ((|+f 2 )=^p/ 2 , or 
(C.+f;)s^g. The spectral domain model of the discontinuity is tlien employed to 
calculate the resonance frequency f„ of the resonator plus discontinuity. The 
terminating stub length f, is changed by a small amount, and then the new length 
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of the stub at the input which is required to re-establish resonance at the same 
frequency f^ as above, is determined. This process is repeated until a sufficient 
number of stub-length pairs ((,,£ 2 ). which produce resonance at the same frequency 
f„. is available for the calculation of the element values of the desired equivalent 
circuit representation. 

In the case when an L-section equivalent circuit model is wanted, it turns out 
that, although only 2 independent circuit elements are included in the equivalent 
circuit, a determination of 3 pairs of (£,,£ 2 ) values is required, in order to effect 
separation of variables in the circuit equations. This is shown in section III. 

Likewise, when an equivalent T- or Pl-section presentation is evaluated for 
an asymmetrical discontinuity, a total of 4 of the (fi.fn) pairs must be resonated. 
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in order to provide the data needed for the determination of the 3 independent 
circuit elements in the equivalent circuit. 
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ni. COMPUTER PROGRAM CONSTRUCTION 


A. CURRENT DISTRIBUTIONS 

In actual computations for the dominant mode, the strip current densities 
and J„, have been chosen to have physically plausible distributions, as shown in 
Fig. 6, for the continuous strip which has width of 2w, and length of 2{. 


1. Continuity Case 
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The coordinate forms of the current distributions in Fig. 6 are: 


^ 2W^ 


(11a) 




(lib) 


1 TT V 

J^(a') = —sin—^ 

H-j U-, 


(11c) 


JA7,)=- 


(lid) 


Fourier transforms of these aie; 


2sin(A-„hi,) 3 

J , {k) -__JLJ_ +_I_ [XIN] 


(12a) 


/ VT 7 T- , 1 . 1 2sin(;c„Kj) 2 [1-cos (ir^h',) ) 

( XIN = cos (kw.) ~ _,_ 4--^_J— ) 

KPl 


2 ] 






47lCOSp^i 




27tsin{/c„Vj) 




cos (p^i) 
Pn-Z' 


sin(p^i) 

(P.i)^ 


And then, forming products of the x- and z- dependent factors, 


J,(A„,p;,)=c7,(;c,j-j^.(p^) 


^^.<^„^P.)=-^(ic„)-J,(p„) 


(12b) 


(12c) 


(12d) 


(13a) 


(13b) 
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2. Step Discontinuity Case 














The coordinate forms of current distributions m Fig. 7 are: 


2W^ Wy 


(14a) 


fl+— 3 

2 (Vj 


(14b) 




(14c) 


T /«\ Z 


(14d) 


J^, (x) =-i-sin— 

Wl U’j 


(14e) 


J .(x) = —sin — 
u’, tv. 


(14f) 


2 if 


(14g) 


2 if 


(14h) 


The Fourier transformed distributions of Eqs. 05) are; 
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^ =Same as Eg. (125) 


(15a) 








[XIN] 


(15b) 


( XIN H CO. ^ 11-cos (;c„>.,) ] 


kh\ 


(^.^z) 


'■^ 21 (p ffl) =Same as Eg. (I2b) 


(15c) 



47CCOS (P„ij) 


(15d) 


XI Eg. (12 c) 


(15e) 




(15f) 
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(16b) 


J_,.. (A., p;,) = c,J 




ik^) -j 




where the individual terms have non-zero values only in the line segment indicated 
by the corresponding subscript, as explained above. 

From Eqs. (5) and (6) 

- - (17a) 

= E E ^ P J ^ zzJ zi ^K’ P 

n=J m=l 
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(17b) 




n=l m'l 


X21‘E E fj',,- p .> z x.j „ P J 


(170) 


n=l w-1 


K22^E Z P .' Z' „j „ (J:,. P „) 


(17d) 


n'2 m=3 
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B. COMPUTATION PROCEDURES 


1. Change of Transform Variable 

In the present work, the strip line resonator was fully enclosed in a metal 
shield. This permitted the use of a finite Fourier transform rather than the integral 
transform which is typical of the infinite-line calculations. 

The inner products can then be carried out as truncated summations, 
saving much computer timefRef.9]. 

p CB 

Jo 

Strictly speaking, the limits are n=-<x5, m=-«> on the low side, unless symmetiy' 
holds, and the sum is: 

The wavenumbers used for the two summations are: 

x-direction 

z-direction 










E 


EE 

/)*1 /P-3 
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As shown on Table II and III, a summation over n of 20 terms is enough, 
and for the summation over m, 100 terms leads to convergence within 1 MHz error 
which is negligible in comparison with the resonant frequency. 


TABLE II VARIATION OF NUMBER OF TERMS N 


n 

m 

"o (GHz) 

10 

100 

5.7168 

20 * 

100 

5.7250 

40 

100 

5.7261 


TABLE III VARIATION OF NUMBER OF TERMS M 


n 

m 

fo (GHz) 

20 

40 

5.7072 

20 

80 

5.7228 

20 

100 ♦ 

5.7250 

20 

200 

5.7275 
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2. Equations for Some Factors 


a. Propagation Constant P 

By use of the program shown in Appendix A, P is computed at each 
resonant frequency and line-width Wj and Wj. This program is constructed using 
the moment method of calculation. 


b. Effective Dielectric Constant 

The calculation of uses P(f), frequency-dependent propagation 
constant, and P„^, quasi-static propagation constant [Ref.3j. 


e = f-0_LLL 12 

eff~ I J 

” air 


(18) 


c. Characteristic Impedance Zff) 

The computation of characteristic impedance uses an empirical 
calculation for the quasistatic value, "Lff)) which has been confirmed to be accurate 
[Ref. 10]. This value is then corrected to the wanted frequency f by use of e^ff in 
the expression :[Ref,4] 




6f f 


(19) 
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for 0 < w < a/2; 


ZJO) = 60[V+J^-ln[ 


1+4 ( 


V = -1.7866 - 0.2035(t/b) + 0.4750(a/b) 
R = 1.0835 + 0.1007(t/b) - 0.09457(a/b) 


for a/2 < w < a; 

Z^(0) = 1207t[V+i?(.^%1.3930 + 0.6670-ln(2l^4l.444) 

b b 

V = -0.6301 - 0.07082(t/b) + 0.2470(a/b) 

R =: 1.9492 + 0.1553(tA)) - 0.5123(a/b) 
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C. PRESENIATION OF THE L-CIRCUIT COMPONENTS WHICH 
REPRESENT THE DISCONTINUITY 

The network is considered as a two-port, i.e., 4-terminal, element tenninated 
in open-ended stub sections, as indicated in Fig. 8. 



Figure 8 L-circuit representation 


The impedance at the input to port 1 is Zj„|. The two stubs have input 
impedances: 


= J ZojCO t P j i j 

(20a) 


(20b) 


The input impedance to the network is then: 


32 







Setting Zjni = z,*, and expanding the expression, we get: 


(22a) 

where we define: A s z,* - Z 2 
Bsz,‘ 

C s - 2 ,* 2 , 

We now assume that we have values A, B, C and A’,B’,C’, as well as 
A",B".C", of the constants defined above, corresponding to different stub-length 
pairs, ((, ,£ 2 ). (£/,£ 20 ) and (£'i",f 2 ^ 0 . all of which have been determined by 
spectral-domain calculations to resonate at the same resonance frequency f„. 

We rewrite the equation above for each of the stated cases: 


Z.rA'-Z^Z^-^Z^B'=C' 

(22b) 

z.^A"-z^z^^z^B"-C" 

(22c) 


Eqs.(22) are solved simultaneously to obtain the values for the circuit 


components z, and 7 .^: 





(23a) 


(C-C") (B-B0-(C-C0 [B-B") 
{A-A") {B-B') - (A-A') {B-B") 


{C-C')-Z^kA-A') 

(B-BO 


(23b) 


The component values for a three-element PI- or T-circuit can be calculated 
in a completely analogous manner, but with use of a set of four ((,,( 2 ) pairs. 
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D. RESULTS OF COMPUTATION 


1. Strip Resonator ■with Uniform Width 

a. Resonant Frequencies with Varying Height of Lower Air Layer 
The lower air gap is varied from the bottom of the shield to the 
middle of the shield height. When the dielectric material is near the bottom, 
comparable to the case of the microstrip line configuration, the resonant frequencies 
of SSSL show the same dependence on width w, as for a microstrip resonator 
(Ref.iJ. As the air gap becomes larger, the tendency is reversed[Fig.9J. 
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b. Resonant Frequencies with Varying Lengths of Stub Segments on 
the Test Resonator 

When the length of strip is greater, the resonant frequency is smaller 
as expected, because the resonant wavelength is proportional to the length of strip. 



Figure 10 Resonant frequency vs. length of strip with chajige of lowei air gap 











c. Effects of Change of Substrate Dielectric Constant 

As the substrate dielectric constant is increased, the resonant 
frequency of the strip resonator is reduced. But the variation of resonant 
frequency with change of lower air gap is faster as dielectric constant 
increases[Fig.ll]. 



Figure 11 Resonant frequency vs. width of strip with change of 
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2. Strip Resonator with Step Discontinuity in Width 


a. Resonant Frequencies with Discontinuity at Center of Resonator 
When the length C|(=C 2 ) is decreased, the resonant frequency is 
increased as in the uniform width resonator shown in Fig. 10. A greater difference 
between widths w, and Wj results in a larger decrease of resonant frequency. 



Figure 12 Resonant frequency vs. Ci(=C 2 ) with different Wj/Wj 
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b. Resonant Frequencies with Varying Ratio of Lengths 


A change of the ratio of the lengths affects the resonant frequencies. 
In turn, the greater difference between the lengths of the segments of different 
widths the higher is the resonator frequency. However, the resonant frequency 
converged to a single value, as the ratio of two lengths is was made larger. 


X 

o 
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Figure 13 Resonant frequency vs. Cj/Cj with different w,/w 2 
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c. Capacitance of Step with varying Parameters 


The equivalent-circuit capacitance of Fig,8 is dependent on the 
resonant frequency. To get the different resonant frequencies the length of stub, 
{i(={ 2 )* is changed. When (2 becomes shorter, there is corresponding C, which 
produces resonance with the same frequency as the reference case. As mentioned 
before, 3 pairs of ((,, Cj) are needed. Here, the reference case,®,=( 2 , can not be 
one of the three because it was found that using this case resulted in negative 
values of capacitance resulting from the computation. More investigation is 
required to determine the reason for this behavior. 

As shown in Fig. 14 & Fig. 15, a greater value of the ratio w,/w 2 
corresponds to a greater value of discontinuity capacitance, as expected[Ref.7]. 
However, when the W'j/Wj approaches to unity, the capacitance does not decrease 
continuously but was found to increase. Since the limiting case, w,/w 2 =l, should 
correspond to zero capacitance, the result obtained seems to be due to the fact that 
the difference terms in Eq.(23) represent the taking of differences between nearly- 
equal quantities, and a small error in the individual quantities will produce a large 
error in the result. Other factors are the approximate nature of the basis functions 
and/or the small number of basis functions. 

Shorter values of f, and ^2 correspond to higher resonant frequency. 
The variation of resonant frequency shows the frequency-dependence of the 
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equivalent- circuit capacitance [Fig. 16]. Also, a greater substrate dielectric constant 
corresponds to a greater equivalent-circuit capacitance as expected on physical 
grounds [Fig. 17]. 
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Figure 14 Capacitance vs. Wj/Wj with change of initial resonator length 



Figure 15 Capacitance vs. Wj/Wj with change of e, 
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d. Inductance of Step with vaiying Parameters 


The procedure for the calculation of the equivalent circuit inductance 
is analogous to that used in the case of the ciq)acitance. The calculated inductance 
did not converge to zero when Wj/Wj goes to unity, also[Ref.8]. The inductance 
does not show a monotonic incmase for Wj/Wj greater than two[Fig. 18]. This 
seems to be due to the approximate nature of the calculation. 

Fig. 19 generally shows similar behavior of the calculated inductance, 
with several different values of dielectric constant. The inductance decreases as the 
resonant frequency increases, however[Fig.20]. Increasing the substrate dielectric 
constant causes a uniform increase in the discontinuity inductance[Fig.21]. 
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Figure 18 Inductance vs. Wj/Wj with change of initial resonator length 



Figure 19 Inductance vs. w/wj with change of e, 
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Figure 20 Inductance vs. resonant frequency with change of w/wj 
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IV. CONCLUSION 


This thesis follows the methods originated by Itoh for the analysis of the 
shielded suspended stripline resonator[Ref.l & 3]. A fully-shielded enclosure was 
employed, allowing the use of a finite Fourier transform in two coordinate 
directions. This change reduced the computation time significantly, while 
producing accuracy comparable with that of the use of the integral transform along 
the line axis. 

The perturbed-iC' onator technique permitted the use of the strip current- 
density distributions suggested by Itoh[Ref.l]. These are known to give accurate 
results in Galerkin’s method. 

We may assume various functional forms for the current distribution on the 
line-stub of the resonator. However, no single distribution can give an exact result, 
though a larger number of distribution functions should lead to more accurate 
results. The choice of distribution function is important because of the very small 
number of functions used in this work. 

The data of Fig.9 shows an interesting change of the slope of the resonant 
frequency vs width ratio with varying height of lower air gap[Fig.l]. Fig. 14 & 18 
shows increasing of capacitance and inductance respectively, as the resonant 
frequency decreases. But, as mentioned above, these graphs show some deviation 
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from expectation, in that the capacitance and inductance values fail to become zero 
with no step in width. This behavior should be given further research. 
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APPENDIX. COMPUTER PROGRAMS USING FORTRAN 


A. RESONANT FREQUENCY 


FILEt STEPFO FOR A1 


INTEGER I,N,M,COMPARE,COUNT 
REALX8 W1,W2,L1,L2,A,D,T,H,ER,F 
REALX8 KN,BETA,OMEGA,RL,RRL 
REALX8 PI,E0,MU0,C 
REALX8 GAM1S,GAM2S 

REALX8 JZIN1,J2IN2,JZKN1,JZBE1,JXKN1,JXBE1,JX1,JZ1 
REALX8 JZKN2,JZBE2,JZBE3,JXKN2,JXBE2 

REALX8 CMP,LIMIT,OLDF,NEHF,LIMF 
C 

LOGICAL LLOG 
C 

COMPL EXX16 GAMAl,GAMA2,GAMAS,CTl,CT2,CT3 
C0MPLEXX16 ZEN,ZENR,ZED,ZHN,ZHD,ZHNR,ZDEN 
C0HPLEXX16 ZE,ZH,ZZZ,ZZX,ZXZ,ZXX 
C0MPLEXX16 K11,K12,K22,J,SK11,SK12,SK22,DET 
CHARACTERX8 FNAME 
CHARACTERX12 SNAME 
C 

0PEN(1,FILE='DATAIN1*,STATUS=»0LD*) 

PRINTX, 'ENTER OUTPUT DATA FILE NAME <FNAME>=' 
READX, FNAME 
0PEN(2,FILE=FNAME) 

SNAME=FNAME//'PLT' 

OPEN(3,FILE=SNAME) 

HRITE(2,x) • FILE NAME = ',FNAME 
WRITE(3,X) ' FILE NAME = »,SNAME 

PARAMETERS 
C=3.E8 

J=CMPLX(0.,1.) 

PI=4.XATAN(1.) 

E0 = 1.E-9/(36 .XPI) 

MU0=4.XPIX1.E-7 
C LIMIT=l.E-2 

LIMIT=l.E-5 
C0MPARE=1 
C0UNT=1 
LL0G=.TRUE. 

C 

CALL READF(W1,W2,L1,L2,A,D,T,H,ER) 

PRINT*,'ENTER THE NUMBER OF OUTER SUM LOOP N = » 
READK,N 
C N=20 

WRITE(2,x) 'THE NUMBER OF OUTER SUM LOOP N =',N 

WRITE(3,X) 'THE NUMBER OF OUTER SUM LOOP N =»,N 

C 

PRINT*, 'ENTER THE NUMBER OF INNER SUM LOOP M = ' 
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READX, M 

C M=20, 40, 100, 500, 1000, 2000, 3000, 4000 ETC. 

WRITE(2,X) 'THE NUMBER OF INNER SUM LOOP M =',M 

WRITE(3,X) 'THE NUMBER OF INNER SUM LOOP M =‘,M 

C 

RL=0.10D0 

RRL=RLX1.D3 

C 

HRITE(2,x) *RL = *,RRL 
WRITE{3,X) »RL = SRRL 
C 

1 PRINTx,'ENTER INITIAL FREQUENCY (GHZ) • 

READX,F 
F=F*1.E9 

PRINTX,'FREQUENCY = ',F 
WRITE(2,X) 'FREQUENCY = •,F 

0MEGA=2.XPIXF 

K11=CMPLX(0.,0.) 

K12=CMPLX(0.,0.) 

K22=CMPLX<0.,0.) 

DO 20 1 = 1,N 
KN=IXPI/A 
KN=(I-.5)XPI/A 

JZIN1=C0S(KNXW1)-2.XSIN(KNXW1)/(KNXW1) 

+ +2.X(1-C0S(KNXW1))/(KNXW1)XX2 

JZIN2=C0S(KNXW2)-2.XS1N(KNXW2)/(KNXW2)+2.X(1-C0S(KNXN2)) 
+ /(KNXW2)XX2 

JZKN1=2.XSIN<KNXH1)/(KNXW1)+3.XJZIN1/(KNXH1)XX2 
JZKN2=2.XSIN(KNXW2)/(KNXW2)+3.XJZIN2/(KNXW2)XX2 
JXKN1=-2.XPIXSIN(KNXW1)/(PIXX2-(KNXW1)XX2) 
JXKN2=-2.XPIXSIN(KNXW2)/(PIXX2-CKNXW2)XX2) 

SK11=CMPLX<0.,0.) 

SK12=CMPLX(0.,0.) 

SK22=CMPLX(0.,0.) 

DO 40 IM=1,M 

BETA=(REAL(IM)-.5)XPI/RL 
BETA=REAL(IM)XPI/RL 

GAM1S=KNXx2+DETAxx2-EOxMU0XOMEGAxx2 

CALL GAMCT(GAMA1,GAM1S,CT1,D) 

CALL GAMCT(GAMA3,GAM1S,CT3,H) 

GAM2S=KNXX2+BETAXX2-EOXMUOXERXOMEGAXX2 

CALL GAMCTCGAMA2,GAM2S,CT2,T) 

ZEN=GAMA2XCT1/ER + GAMA1XCT2 
ZENR=-J/(OMEGAXE0) 

ZED=CT1XCT2+CT1XCT3XGAMA2/(GAMA3XER) 

+ +CT2XCT3 + ERXGAMA1/GAMA2 

ZE=ZENRXZEN/'ZED 
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ZHN=GAMA1XCT1 +GAMA2XCT2 
ZHNR=JXOMEGAXMUO 

ZHD=GAMA1XCT1KGAMA2XCT2 + GAMA1KCT1XGAMA3XCT3 
+ + GAMA2XCT2XGAMA3XCT3 + GAMA2XX2 

ZH=ZHNRXZHN/ZHD 
C 

ZDEN=KNXX2 + BETAXX2 
ZZZ=-(ZExBETAXK2 + ZHXKNXX2)/ZDEN 
ZZX=-KNXBETAX(ZE-ZH)/ZDEN 
ZXZ=ZZX 

ZXX=-(ZEXKNXX2 + ZHXBETAXX2)/ZDEN 
C 

JZBE1=A.xPIxC0S(BETAxL1)/(PIXx2-(2.xBETAxL1)xx2) 

JZBE2=A.XPIXC0S(BETAXL2)/(PIXX2-(2.XBETAXL2)XX2) 

J21=JZKN1XJ2BE1+JZKN2XJZBE2 

JZ1=JZ1X0.5D0 

JXBEl=C0S<BETAxLl/2.)/{BETAxLl/2.) 

+ - SIN(BETAXLl/2.)/(BETAxLl/2.)xx2 

JXBE2=C0S(BETAXL2/2.)/(BETAXL2/2.) 

+ - SIN<BETAXL2/2.)/(BETAxL2/2.)xx2 

JX1=JXKN1XJXBE1+JXKN2XJXBE2 

C 

SK11=SK11+JZ1XZZZXJZ1 
SK12=SK12+JZ1XZZXXJX1 
SK22=SK22+JX1XZXXXJX1 
AO CONTINUE 

WRITE(2,X) JZ1,JX1 

K11=K11+SK11 
K12=K12+SK12 
K22=K22+SK22 
CONTINUE 

PRINTX, 'CHECK Kll = ',K11 
WRITE(2,X) 'CHECK Kll = •,K11 
PRINTX, 'CHECK K12 = ',K12 
WRITE(2,X1 'CHECK K12 = ’,K12 
PRINTX, 'CHECK K22 = ',K22 
WRITE(2,X) 'CHECK K22 = ’,K22 


DET=K11XK22-K12XK12 
PRINTX, 'DET = '/DET 
HRITE( 2 ,xi 'DETERMINENT = ',DET 


CALL FDS0RT(F,DET,COUNT) 
C0UNT=C0UNT+1 


IF (LLOG) THEN 

PRINTX, 'MAKE SURE OPPOSITE 
PRINTX, 'IF SAME SIGNS THEN 


SIGNS BETNEEN 1ST 8 2ND DETS' 
ENTER 9 , OR ENTER O' 
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READX, UK 
IF (IJK.EQ.9) THEN 
0LDF=F 
GOTO 1 
ENDIF 

LL0G=.FALSE. 

ENDIF 


CMP=REAL(DET) 

LIMF»ABS(F-OLDF) 

PRINTX, 'CHECK LIMIT FREQ »,LIMF 
WRITE(2,K) 'CHECK LIMIT FREQ = '^IMF 

IF (ABSCCMP).GT.LIMIT.AND.LIMF.GT.1.E3) THEN 
IF (ABS<CMP).GT.LIMIT.AND.LIMF.GT.l.E-1) THEN 
CALL SLVAL(CMP,OLDF,F/NEWF,COMPARE) 

OLDF=F 
F=NEWF 
GOTO 11 
ENDIF 

PRINT*, 'LAST FREQUENCY = ',F 
WRITE(2,K) 'LAST FREQUENCY = ',F 

CALL FDSORT<F,DET,100) 

STOP 

END 


C****XX***K*X*****X)(XXXXXX**XXX*XXXX*XXXXXX»(*X*XXXXXK*XXXX**X 
SUBROUTINE READF(W1,W2,L1,L2,A,D,T,H,ER) 

REALX8 Wl,W2,Ll,L2,A,D,T,H,ER 
INTEGER DDDD 

PRINTx,'DATA INPUT i KEYBORD(I),FILE(0)' 

READX,DDDD 

9 IF(DDDD.EQ.O) THEN 

READ(1,X) W1,W2,L1,L2,A,D,T,H,ER 
ELSEIF(DDDD.EQ.l) THEN 

PRINT*,'ENTER THE HI 8 W2' 

READ*,HI,H2 

PRINT*,'ENTER LI & L2' 

READ*,LI,L2 
PRINT*,'ENTER A' 

READ*,A 

PRINT*,'ENTER D, T, 8 H» 

READ*,D,T,H 
PRINT*,'ENTER ER' 

READ*,ER 

ELSE 

PRINT*,'SELECT I OR 0 ' 

READ*, DDDD 
GO TO 9 
ENDIF 

PRINT 32, W1,W2,L1,L2,A,D,T,H,ER 
WRITE(2,32) W1,W2,L1,L2,A,D,T,H,ER 
WRITE(3,32) W1,W2,L1,L2,A,D,T,H,ER 
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W1 » WlKl.E-3/2. 

H2 = W2X1.E-3/2. 

LI = LlXl.E-3 
L2 = L2X1.E-3 
A = AKl.E-3/2. 

D = Dxl.E-3 
T = Txl.E-3 
H = HKl.E-3 

32 F0RMAT(T3,'UNITtMILIMETER' 

+ ,/,T3,'Wl = •,F9.5,5X,•M2 = •,F9.5 

+ ,/,T3,'Ll = •,F9.5,5X,*L2 = »,F9.5 

+ ,/,T3,'A = •,F9.5 

+ ,/,T3,'D = •,F9.5,5X,'T = •,F9.5,5X,'H = ',F9.5 

+ ,/,T3,'ER = •,F9.5,/) 

RETURN 

END 

CXXXXXKXXXKXXKXXXXiOtXKXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXKK 

SUBROUTINE GAMCTCGAMA,GAMS,CT,LENGTH) 

C0MPLEXX16 GAMA,CT 

REALX8 GAMS,LENGTH,GAMR,CTR 

GAMR =DSQRT(DABS(GAMS)) 

IF(GAMS.LT.O.DO) THEN 

GAMA = CMPLX(0.D0,GAMR) 

CTR = -l./TAN(GAMRXLENGTH) 

CT = CMPLX(0.D0,CTR) 

ELSE 

GAMA = CMPLX(GAMR,0.D0) 

CTR = l./TANHCGAMRXLENGTH) 

CT = CHPLX(CTR,0.D0) 

ENDIF 

RETURN 

END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE SLVAL(FUN,OXV,XV,NEWXV,SN) 

INTEGER SN 

REALX8 FUN,OXV,XV,NEMXV,POSXV,NEGXV 

IF(SH.EQ.l) THEN 
IFCFUN.GT.O.) THEN 
POSXV = XV 
NEGXV = OXV 
ELSE 

NEGXV = XV 
POSXV = OXV 
ENDIF 

NEWXV = (POSXV + NEGXV)/2. 

SN = SN +2 
ELSE 

1F(FUN.GT,0.) THEN 
POSXV = XV 
ELSE 

NEGXV = XV 
ENDIF 

NEWXV = (POSXV +NEGXV)/2. 

ENDIF 

RETURN 

END 



CXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXKXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXX' 

SUBROUTINE FDSORK FQY,SDET,ICOUNT) 

INTEGER ICOUNT/IC 

REALXB FREQY(100),FQY,RTEMP,MAGDET,TFQY 
C0MPLEXX16 DETMT(100),SDET,CTEMP,TSDET 
TFQY = FQY 
TSDET = SDET 
IF(ICOUNT.EQ.l) THEN 
FREQY(l) = TFQY 
DETMT(l) = TSDET 
ELSEIF(ICOUNT.NE.IOO) THEN 
DO 10 I=1,1C0UNT-1 

RTEMP = FREQY(I) 

CTEMP = DETMT(I) 

IF (RTEMP.GT.TFQY) THEN 
FREQY(I) = TFQY 
TFQY = RTEMP 
DETMTd) = TSDET 
TSDET = CTEMP 
ENDIF 

10 CONriNUE 

FREQY(ICOUNT) = TFQY 
DETMTdCOUNT) = TSDET 
IC = ICOUNT 

ELSEIF(ICOUNT.EQ.IOO) THEN 
WRITE(3,110) 

WRITE(3,120) 

no F0RMAT(//,T8,'FREQUENCY ST23,'MAG OF DET ', TA3,'DETERMINENT') 
120 FORMAT<1X,TAO,'RE.•,T56,'IM.*) 

DO 20 J =1,IC 

MAGDET *DSQRT(REAL(DETMT<J))XX2+DIMAG(DETMT(J))XX2) 

HRITEC 3,210) FREQYCJ),MAGDET,DETMTCJ) 

210 F0RMAT<T5,F15.2,T20,F15.3,T58,2(E15.8)) 

20 CONTINUE 

ENDIF 
RETURN 
END 
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B. PROPAGATION CONSTANT P 


FILE. BETA FOR A1 


INTEGER MMAX,NMAX 

REALX8 N.A,D,T,H,PI/ER/F6,OMG,EO,KN 

REALMS G1S2,B0,B00,DBTA,FLAG»AA2,CC,WN 

REALXB G1,CT01,CT03,GT1,GT3,G2S2,G2,CT02,KE1,KH1,J21 

REALK8 B,OUT,EREFF,FC,BIM,Q,R 

DIMENSION B(60),OUT(30) 

C0MPLEXX16 GM1,GM3,CT1,CT3,GM2,CT2,DENE,ZE,KE,DENH,ZH,KH 
COMPLEX TRMN,SUMUP/JZ11 
C 

0PEN(2,FILE=*BETAVAL*,STATUS=‘0LD») 

C 

PI=4.D0XATAN(1.D0) 

FLAG=0 

C 

MMAX=30 

NMAX=30 

C 

PRINTX,‘ENTER FREQ.,GHZ.' 

READJ<, FG 
C 

PRINTX,‘ENTER SHLD. WDTH(MM>' 

REAO)(, AA2 
A=AA2/2000 

PRINTX,'ENTER SHLD. HT.(MM)' 

READX,CC 
CC=CC/1.D3 

PRINT*,‘ENTER LOWER AIRGAP,(MM)‘ 

READ*,D 

D=D/1.D3 

PRINT*,'ENTER SUBSTRATE THICKNESS(MM)' 

READ*,T 

T=T/1.D3 

H=CC-(D+T) 

PRINT*,'ENTER THE SUBSTRATE DIELEC. CONST. ' 

READX,ER 

PRINT*,'ENTER THE LINEWIDTH(MM)» 

READ*,WW 

W=WW/2.D3 

Q=0.2077+1.2177XT/CC-0.0836AKAX2./CC 
PRINT*,'0=',Q 

R=0.03451-0.1031*T/CC+0.01742XA*2./CC 

PRINT*,'R=',R 

TT=Q-RXL0G(W*2./CC) 

SS=LOG(l./SQRT(ER)) 

EREFF=1./(1+TT*SS) 

EREFF=EREFF**2 
PRINT*, 'EREFF = ',EREFF 


0MG=FGXPIX2E9 
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E0=8.8<iD-12 

C BO«(OMG/3D8)KSQRT(<ER+l)/2) 

BO=(OMG/3D8)KSQRT(EREFF) 

BOO-O.2OD0 
DBTA^O.IODO 
GO TO 4 
5 BOBBIN 

BOO=0.85DO 
DBTA’^O.OOTDO 
A DO 2 M=1,MMAX 
SUMUP=0 

B(M)-BOX(BOO-t-MxDBTA) 

DO 1 N=1,NMAX 
KN=CN-0.5D0)KPI/A 
G1S2=KNXX2+B(M)XX2-0MGXK2/9D16 
G1=SQRT(ABS(G1S2)) 

IF(GlS2.LT.O) THEN 
GM1=CMPLX(0.D0>G1) 

GM3SCM1 

CT01*-1.DO/TAN(01XH) 

CT03=-1.D0/TAN(G1XD) 

CT1=CMPLX(0.D0,CT01) 

CT3=CMPLX<0.D0,CT03) 

ELSE 

GM1=CMPLX(G1,0.D0) 

GT1=1.D0/TANH(G1XH) 

GM3=GM1 

GT3=1.D0/TANH(G1XD) 

CTl=CMPLX<GTl,O.DO) 

CT3=CMPLX(GT3,0.D0) 

ENDIF 

C 

G2S2=KNXX2+B<M)XX2-ERX0MGXX2/9.D16 
G2=SQRT(ABS(G2S2)) 

IF (G2S2.lt.0) THEN 
GM2=CMPLX(0.D0,G2) 

CT02=-1.D0/TAN<G2XT) 

CT2=CMPLX(0.D0,CT02) 

ELSE 

GM2=CMPLX(G2,0.D0) 

GT2=1.D0/TANH(G2XT) 

CT2=CMPLX(GT2,0.) 

ENDIF 

C 

DENE=CT2XCT3+CT1XCT3X(GM2/GM1)/ER+CT1XCT2+{GM3/GM2)XER 

ZE=(GM2XCT3/ER+GM3XCT2)/DENE 

KE1=-1.DO/(OMGXEO) 

KE=CMPLXCO.DO,KEI) 

ZE=KEXZE 

DENH=GM1XCTIXGM2XCT2+GM1XCTIXGM3XCT3+GM2XCT2XGM3XCT3+GM2XX2 

ZH=(GM2XCT2+GM3XCT3)/DENH 

KH1=0MGXPIXAD-7 

KH=CMPLX(O.DO,KHI) 

ZH=ZHXKH 

JZ1=COSCKNXW)-2.DOXSIN(KNXW)/(KNXW)+2.DOX(1.-COSCKNXW))/(KNXW)XX2 

JZI=2.D0XSIN<KNXH)/(KNXW)+3.D0/((KNXW)XX2)XJZ1 

JZII=CMPLX<JZ1,O.DO) 

TRMN=-(B(M)XX2XZE+KNXX2XZH)/(B(M)XX2+KNXX2)J;JZ11XX2 
SUMUP=TRMN+SUMUP 
1 CONTINUE 
C 
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OUT<M)sCABS(SUMUP) 

CONTINUE 


DO 3 M=2,(MMAX-1) 

IF <(OUT(M).LT.OUT(M-l)).AND.(OUT(M).LT.OUTCM+l))) THEN 
BIN=B(M) 

ELSE 

ENDIF 

CONTINUE 

IF CFLAG.EQ.O) THEN 
FLAG=1.D0 
GO TO 5 
ELSE 

PRINTX,'SUSP. STRIPLINE BETA = SBIN 
WRITE(2,X) ‘SUSP. STRIPLINE BETA = •,BIN 
ENDIF 
STOP 


oo 


C. EQUIVALENT CIRCUIT CAPACITANCE & INDUCTANCE 


FILE. EQICIR FOR A1 


INTEGER I 

REALX8 BETA1,BETA2»Z01,Z02 
REALX8 PI,F,OMEGA,Wl,M2 

COMPL EXX16 A, B, C, J, Z3,24, NUMZ5, ;,ENZ3, NUMZ4, DENZ4 
C0MPLEXX16 CAPAC,INDUC,MC,ML 
DIMENSION A(3),B(3),C(3) 

CHARACTERX8 FNAME 

PRINTX, 'ENTER OUTPUT DATA FILE NAME' 

READX,FNAME 
0PEN(2,FILE=FNAME) 

J=CMPLX<0.D0,1.D0) 

PI=4.D0XATAN(1.D0) 

PRINTX,'ENTER RES. FREQUENCY (GHZ)' 

READX,F 

WRITE(2,X) 'RESONANT FREQUENCY(GH2) =',F 

F=FX1.E9 

OMEGA=2.DOXPIXF 

PRINTX,'ENTER Ml 8 M2' 

REA0X,M1,M2 
MRITE(2,130) Ml,M2 

130 F0RHAT(T3,'M1 = ',F3.1,'MM*,5X,'M2 = ',F8.4,'MM') 
PRINTX, 'ENTER ZOl & Z02' 

READX, 201,202 
PRINTX,'ZOl = •,Z01 
PRINTX,*Z02 = ',Z02 
MRITE<2,100) Z01,Z02 

100 FORHAT(T3,'Z01 = ',F8.4,5X,»Z02 = »,F8.4) 

PRINTX, 'ENTER BETAl 8 BETA2’ 

READX, BETAl,BETA2 
PRINTX, 'BETAl = ',BETAl 
PRINTX, 'BETA2 = ',BETA2 
MRITE(2,110) BETAl,BETA2 

no F0RMAT(T3,'BETAl = ',F8.4,5X,'BETA2 = ',F8.4) 

DO 10 1=1,3,1 

CALL Z1Z2ABC(BETAI,BETA2,Z01,Z02,J,I,A,B,C) 

10 CONTINUE 

NUMZ5=(C(1)-C(3))X(BC1)-B(2))-(C(1)-C(2))XCB(1)-B(3)) 
DENZ3=<A(1)-A(3))X(B(1)-B(2))-(A(1)-A(2))X(B(1)-B(3)) 
NUMZ3=(C(1)-C(2))X(B(2)-B(3))-(C(2)-C(5))X(B(1)-B(2)) 
DENZ3=(A(1)-A(2))X(B<2)-B(3))-CA(2)-A(3))X(B(1)-B(2)) 
Z3=NUMZ3/DENZ3 

NUMZ4=(C(1)-C(2))-Z3X(A(1)-A(E)) 

DENZ4=B(1)-B(2) 

Z4=NUMZ4/DENZ4 

PRINTX,'NUMZ3 =',NUMZ3,'DENZ3 =',DENZ3 
PRINTX,»NUMZ4 =',NUMZ4,*DENZ4 =',DENZ4 
MRITE(2,X) 'NUMZ3 =',NUMZ3,'DENZ3 =',DENZ3 
MRITE(2,X) 'NUMZ4 =',NUMZ4,'DENZ4 =',DENZ4 
PRINTX,»Z3 = ',Z3 
PRINTX,*Z4 = ',Z4 
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WRITE(2,X) »23 = ',23 
WRITE(2,X) = ',24 
CAPAC=l./f23XJX0ME6A) 
INDUC=24/(JX0MEGA) 

PRINTX, 'CAPAC.s •,CAPAC 
PRINTX, 'INDUC.s '.INDUC 
WRITE(2,X) 'CAPAC. = •,CAPAC 
WRITE(2,X) 'INDUC. = *,INDUC 
WC=1.DO/(OMEGAXCAPAC) 
NL=OMEGAXINDUC 
PRINTX, 'l/WC = •,WC 
PRINTX, 'WL = ',WL 
WRITE(2,X) »1/WC = »,HC 
WRITE(2,X) 'WL = *,WL 
STOP 
END 


SUBROUTINE 2122ABC(BETA1,BETA2,201,202,J,I,A,B,C) 
INTEGER I 

REALX8 BETA1,DETA2,201,202,L1,L2 
COMPLEXXl^ A,B,C,J,21,22 
DIMENSION A<3),B<3),C(3) 


PRINTX, 'ENTER LI 8 L2 (MM)' 

READX, L1,L2 
PRINTX, 'LKMM) = ',L1 
PRINTX, »L2(MM) = ',L2 
HRITE(2,120) L1,L2 

120 F0RMAT(T3,'LKMM) = », F8.4,5X, • L2(MM) = ',F8.4) 
L1=L1/1.D3 
L2=L2/1.D3 

21=JX201/TAN<BETA1XL1) 

Z2=-JX202/TAN(BETA2XL2) 

PRINTX,»21 = ',21 
PRINTX,'22 = ',22 
WRITE(2,X)'21 = ',21 
WRITE(2,X) '22 = ',22 
A(I)=C0NJG(21)-22 
B(I)=C0NJG<Z1) 

C(I)=-C0NJG(Z1)X 22 

PRINTX,'A(',1,' ) = »,A(I) 

PRINTX,'B(',1,' ) = »,B(I) 

PRINTX,'C<',1,' ) = ',C(I) 

HRITE(2,X) »A(',I,» ) = »,A(I) 

WRITE(2,X) 'BC,!,' ) = ’,B(I) 

HRITE(2,X) 'C(',I,' ) = *,C(I) 

RETURN 

END 
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